Astronomy & Astrophysics manuscript no. G1176 


© ESQ 2008 


September 4, 2008 





The HARPS search for southern extra-solar planets ^ 

XIV. Gl 176b, a super-Earth rather than a Neptune, and at a different period 

T. Forveille^ X. Bonfils^-''^ X. Delfosse^ M. 01110114, 3 udry^, F. Bouchy^ C. Lovls^, M. Mayor^, F. Pepe'*, C. 

Perrier\ D. Queloz"^, N. Santos^, and J.-L. Bertaux^ 



00 

o 
o 

C>0 



o 



> 

o 
in 

q 

On 
O 
00 

o 

> 

• ^ 

S3 



' Laboratoire d'Astrophysique de Grenoble, Observatoire de Grenoble, Universite Joseph Fourier, CNRS, UMR 571 Grenoble, 
France 

- Centro de Astroffsica, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 

^ Centro de Astronomia e Astroffsica da Universidade de Lisboa, Observatorio Astronomico de Lisboa, Tapada da Ajuda, 1349-018 
Lisboa, Portugal 

* Observatoire de Geneve, Universite de Geneve, 51 ch. des Maillettes, 1290 Sauvemy, Switzerland 
' Institut d'Astrophysique de Paris, CNRS, Universite Pierre et Marie Curie, 98bis Bd Arago, 75014 Paris, France 
Service d'Aeronomie du CNRS, BP 3, 91371 Verrieres-le-Buisson, France 

ABSTRACT 

Context. A 10.24 days Neptune-mass planet was recently announced to orbit the nearby M2 dwarf Gl 176, based on 28 radial velocities 

measured with the HRS spectrograph on the Hobby-Heberly Telescope. 

Aims. 

Methods. We obtained 57 radial velocities of Gl 176 with the ESO 3.6m telescope and the HARPS spectrograph, which is known for 
its sub-m/s stability. The median photon-noise standard error of our measurements is 1 . 1 m/s, significantly lower than the 4.7 m/s of the 
HET velocities, and the 4 years period over which they were obtained has much overlap with the epochs of the HET measurements. 
Results. The HARPS measurements show no evidence for a signal at the period of the putative HET planet, suggesting that its 
detection was spurious. We do find, on the other hand, strong evidence for a lower mass 8.4 M^o,,;, planet, in a quasi-circular orbit 
and at the different period of 8.78 days. The host star has moderate magnetic activity and rotates on a 39-days period, which we 
confirm through modulation of both contemporaneous photometry and chromospheric indices. We detect that period as well in the 
radial velocities, but it is well removed from the orbital period and no cause for confusion. 

Conclusions. This new detection of a super-Earth (2 M£„„/,< M sin(i) < 10 M^n,,;,) around an M dwarf adds to the growing evidence 
that such planets are common around very low mass stars: a third of the 20 known planets with M sin(i) < 0.1 Mj„,,and 3 of the 7 
known planets with M sin(i) < 10 M^orrtOrbit an M dwarf, in contrast to just 4 of the ~300 known Jupiter-mass planets. 

Key words. Stars: individual: Gl 176 - Stars: planetary systems - Stars: late-type - Techniques: radial-velocity 



1. Introduction 

Of the ~250 planetary systems currently known from radial 
velocity monitoring, just half a dozen are centered around M 
dwarfs (M<O.6M0) This in part reflects a selection bias, 
since an order of magnitude fewer faint M dwarfs are searched 
for planets than brighter solar-type stars are, but M dwarfs 
also seem to genuinely have fewer massive planets (~ M j„p) 
than the more massive solar- type stars do ( Bonfils et al.||2006 



Johnson et al.|2007 1. They seem on the other hand ( Bonfils et al. 



2006) to have larger numbers of the harder to detect Neptune 
mass and super-Earth planets: a third of the ~20 planets with 
M sin(i) < 0. lMy„pknown to date orbit an M dwarf, in spite of 
solar-type stars outnumbering those by an order of magnitude in 
planet search samples. As a consequence of their small overall 
number, each individual M-dwarf planetary system still plays a 
significant role in defining these emerging statistical properties. 

Send offprint requests to: T. Forveille, e-mail: 
Thierry . Forveille . uj f -grenoble . fr 

* Based on observations made with the HARPS instrument on the 
ESO 3.6-m telescope at La Silla Observatory under programme ID 
072.C-0488 

' http://exoplanet.eu/catalog-RV.php 



Very recently, 'Endl et al. ( 2008 1 announced the discovery 
of a planet with a minimum mass of M sin(i) - 25 M^or//, 
in a 10.24 days orbit around a nearby M2.5 dwarf. Gl 176 
(also HD 285968, HIP 21932, LHS 196) is a V=9.97 
( Upgren,,1974,) member of the immediate solar neighborhood 
(par= 106.2+2.5 mas, d = 9.4 pc. Ferryman & ESA |1997l ). The 
2MASS photometry (Skruts kieTt al.|2006| ) and the parallax re- 
sult in an absolute magnitude of M/^j = 5.74, and together with 
the K-band mass-luminosity relation of Delfosse et al. (2000 ) in 
a mass of 0.50 Mq. Based on the |Bonfils et al. (2005 ) photomet- 
ric metallicity calibration, [Fe/H] is -0.1+0.2 and therefore solar 
within its uncertainty. 

We have independently been monitoiing the radial velocity 
of Gl 176 using the HARPS spectrograph on flie ESO 3.6-m 
telescope, over a period which largely overlaps the epochs of 
the|Endl et al.|([2008) observations. Section 2. describes those 
independent measurements and concludes that they do not con- 
firm the 10.24 days planet. Section 3. takes a closer look at those 
mesurements and finds that they contain two coherent signals, 
with periods of 8.78 and 40.0 days. Section 4 discusses differen- 
tial photometry and variation of chromospheric indices, to con- 
clude that the 40 days signal reflects the stellar rotation period. 
The 8.78 days period on the other hand is due to a bona-fide 
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Table 1. Observed and inferred stellar parameters for Gl 176 



Parameter 




Gl 176 


Spectral Type 


M2V 


V 




9.97 ± 0.03 


n 


[mas] 


106.16 ±2.51 


Distance 


[pc] 


9.42 ± 0.22 


My 




10.10 ±0.06 


K 




5.607 ± 0.034 


Mk 




5.74 ±0.06 


U 


[Lo] 


0.022 






3.5.10-^ 


V sin i 


[kms-i] 


<0.8 


[Fe/H] 




-0.1 ±0.2 




[Mo] 


0.50 



planet, with a minimum mass of only 8.4 M^a,.,/,. Section 5 con- 
cludes with a brief discussion of the new planet. 

2. HARPS Doppler measurements and orbital 
analysis 

We observed Gl 176 with HARPS (High Accuracy Radial ve- 
locity Planet Searcher) as part of the guaranteed-time pro- 
gram of the instrument consortium. HARPS is a high-resolution 
(R - 115 000) fiber-fed echelle spectrograph, optimised for 
planet search programmes and asteroseismology. It is the most 
precise spectro-velocimeter to d ate, with a long-ter m instrumen- 
tal RV accuracy under 1 ms~' (jJVIayor et al. 2003; Santos et al. 



2004 [ Lovis et al. 2005). For ultimate radial velocity precision 



HARPS uses simultaneous exposures of a thorium lamp through 
a calibration fiber. When observing M dwarfs however, we rely 
instead on its very high instrumental stability (nightly instru- 
mental drifts < 1 m s"'). Most M dwarfs are too faint for us 
to reach the stability limit of HARPS within realistic integra- 
tion times, and dispensing with the simultaneous thorium light 
produces much cleaner stellar spectra, suitable for quantitative 
spectroscopic analyses. 

For the V = 9.97 Gl 176 we use 15 mn exposures, and the 
median S/N ratio of our 57 spectra is 60 per pixel at 550 nm. The 
radial velocities (Table|3] only available electronically) were ob- 
tained with the standard HARPS reduction pipeline, based on 
cross-correlation with a stellar mask and a precise nightly wave- 
length calibration from ThAr spectra ( jLovis & Pepe 2007 ). They 
have a median internal error of only 1.1 ms"', which includes 
both the nightly zero-point calibration uncertainty (~ 0.5 m s ') 
and the photon noise, computed from the full Doppler informa- 
tion content of the spectra (Bou chy et al.|200T 



The computed velocities exhibit an rms dispersion of 
5.3 ms This is much above the 1 ms ' internal errors and 
significantly more than we observe for stars with similar chro- 
mospheric activity, but less than the ~8 ms~' expected f rom 



Endl et al. 



(2008) or- 



the 11.7 ms"' velocity amplitude of the 
bit. Fig [T] confirms that the HARPS velocities are more tightly 
packed than both the HET measurements (top panel) and the pre- 
dictions of the lEndl et al. (2008) orbit (lower panel). Its lower 



panel demonstrates that they do not phase on the |Endl et aL 
( |2008| l period, and we verified that the subset of the HARPS 
dataset which overlaps the published HET measurements does 
not either. Since any instrumental or astrophysical noise can only 
increase the velocity dispersion, never decrease it, the HARPS 
measurements set a ~7.5 ms"' ceiling on the radial velocity 
amplitude of a Keplerian orbit (except for unrealistically high 
eccentricities). This forces us to conclude that the |Endl et aL| 
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Fig. 1. Top panel: HARPS (red filled symbols) and |Endl et al. 
(2008) (blue empty symbols) radial velocities of Gl 176 as a 



function of time, overlaid with the prediction of the |Endl et aL 
( |2008| l orbit. Bottom panel: HARPS radial velocities phased at 
the 10.24 days period of the Endl et al. (2008i orbit, overlaid 
with the radial velocity prediction for that orbit. 




10 100 
Period [day] 

Fig. 2. Lomb-Scargle periodogram of the raw HARPS radial 
velocities (top panel), and of the velocities after subtraction of 
the 40-day signal (bottom panel). 



(2008 1 orbit must be spurious, though we do not have a ready 
explanation for why. 



Forveille et al.: A super-Earth planet around the nearby M dwarf Gl 176 



3 



Table 2. Orbital elements for the two-keplerian orbital model of 
Gl 176. 



Element 


Value 


Standard error 


1 [m s 


r 


26.4105 km/s 


0.0004 


Pi [days] 


8.7836 


0.0054 


anet 


ei 


0.0 


Fixed 


"a. 


omi [deg.] 


0.0 


Fixed 




TOi Udb] 


2454399.79 


0.33 


Kli [nVs] 


4.12 


0.52 




Pi [days] 


40.00 


0.11 






0.0 


Fixed 




om2 [deg.] 


0.0 


Fixed 




TO, [jdb] 


2454291.07 


1.31 


1 

CO 


Kl, [m/s] 


4.23 


0.53 





3. Orbital analysis 

Our radial velocity measurements do show coherent structure, 
and a Lomb-Scargle periodogram (Press et al. 1992) shows two 
narrow peaks around 8.8 and 40 days (Fig.|2J top panel). The two 
peaks have similar false alarm probabilities of 0.1%, and their 
spacing is well removed from any significant feature in the win- 
dow function. We therefore analysed them simultaneously and 
searched for 2-planet Keplerian solutions with Stakanof (Tamnz, 
in prep.), a program which uses genetic algorithms to effi- 
ciently explore the large parameter space of multi-planet mod- 
els. Stakanof robustly converged to a 2-keplerian solution with 
periods which match the two periodogram peaks. Subtracting 
the longer period signal from the velocities increases the signif- 
icance of the 8.8 day period in the periodogram (Fig.|2] lower 
panel), further increasing our confidence that this signal is real. 
Subtracting the short period signal, on the other hand, produces 
a periodgram (not shown) with a less convincing 40 days peak. 
The 2-planet model describes our measurements well, but 



certainly not perfectly (cr = 2.5 ms 



2.46 per degree 



of freedom). A Lomb-Scargle periodogram of the residuals of 
this 2-planet solution however shows no significant peak, and the 
significant residuals therefore contain no immediate evidence for 
an additional component. 

Both Keplerian signals have amplitudes of ~4 ms \ which 



with hindsight is well under the sensitivity limit of Endl et aL] 
( 2008| l. Neither of their eccentricities is significant, and we 
therefore adopt circular orbits as our preferred solution (Table 2] 
Fig. [3]i; that choice does not affect any of our conclusions. The 
inner and outer planets, in a Keplerian interpretation of the ra- 
dial velocity variations, have minimum masses (m sin /) of 8 and 
14 MEanh^nd projected semi-major axes of 0.066 and 0.18 AU. 



4. Activity analysis 

Apparent Doppler shifts unfortunately do not always originate 
in the gravitational pull of a companion, because stellar surface 
inhomogeneities, such as plages and spots, can break the bal- 
ance between light emitted in the red-shifted and the blue-shifted 
parts of a rotating star These inhomogeneities then translate into 
rotationally modulated changes of both the shape and the cen- 
troid of spectral lines (e.g. Saar & Donahue|1997[|Queloz et al.| 
[2001 h The activity level of Gl 176 is similar to that of Gl 674 
(Fig. 4jl, where a spot is re sponsible for a 5 ms ' radial velocity 
ignal (|B 
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Fig. 3. Top two panels: Radial velocity measurements phased to 
each of the two periods, after subtraction of the other compo- 
nent of our best 2-planet model. Third panel: Residuals of the 
best 2-planet fit as a function of time (O-C, Observed minus 
Computed). Bottom panel: Lomb-Scargle periodogram of these 
residuals. 



For well resolved rotational broadenings, correlated varia- 
tion in the shape, parametrized by the line bisector, and in the 
centroid, provide an excellent diagnostic of such apparent veloc- 
ity variations. We however measure from our Gl 176 spectra a 
rotational velocity of v sin / < 0.8 km s"' . This small rotation ve- 
locity removes much of the usual power of the bisector test, since 
the bissector span scales with a much higher power of v sin / than 
the centroid ( jSaar & Donahue} 1 997,, Bonfils et al.|2007| l. 

Spots fortunately also produce flux variations, and they typ- 
ically impact spectral indices, whether designed to probe the 
chromosphere (to which photospheric spots have strong mag- 
netic connections), or the photosphere (because spots have 
cooler spectra). We therefore investigated the magnetic activity 
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Fig. 4. Emission reversal in the Ca ii H line in the average spec- 
tra of Gl 674 (M2V, top), Gl 176 (M2.5V, middle), and Gl 581 
(M3V, bottom). Within our 100 M dwarfs sample, Gl 581 has 
one of the weakest Ca ii emission and illustrates a very quiet M 
dwarf Gl 674 and Gl 176 have much stronger emission and are 
both moderately active. 
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Fig. 5. Upper panel: Differential photometry of Gl 176 as a func- 
tion of time. The star clearly varies on a 40-50 days time scale 
with a ~1.3% peak to peak amplitude 



of Gl 176 through photometric observations ( ^4.1| l and detailed 
examination of the chromospheric features in the clean HARPS 



spectra (§4.2 1. 



4.1. Photometric variability 

We obtained photometric measurements with the EulerCAM 
CCD camera of the Euler Telescope (La Silla) during 21 nights 
between November lO"" 2007 and January ll"' 2008 October 
Gl 176 was observed through an Ic filter, to maximize the flux 
of both Gl 176 and a M star in the 11.7' field of view which 
we planned to use as photometric reference. That planned ref- 
erence however proved variable, and we had to fall back to the 
average of two fainter blue stars, with a summed flux of only 
7% of that of Gl 176. In retrospect, this filter choice was there- 
fore suboptimal. To minimize atmospheric scintillation noise we 
took advantage of the low stellar density to defocus the images to 
FWHM ~ 8", so that we could use longer exposure times. The 
increased read-out and sky background noises from the larger 
synthetic aperture which we had to use remain negligible com- 
pared to both stellar photon noise and scintillation. 

We gathered 5 to 7 images per night with a median expo- 
sure time of 3 1 seconds, except on December 29''' when we ob- 
tained sets of 5 images at three well spaced airmasses to measure 



the differential extinction coefficient. We tuned the parameters 
of the Iraf Daophot package ( Stetson|[T987 1 and optimised the 
set of reference stars to minimise the average dispersion in the 
Gl 176 photometry within the individual nights. These param- 
eters were then fixed for the analysis of the full data set. The 
nightly light curves for Gl 176 were normalized by that of the 
sum of the references, clipped at 3-cr to remove a small num- 
ber of outliers, and averaged to one measurement per night to 
examine the long term photometric variability of Gl 176. Gl 176 
clearly varies with a ~1.3% peak-to-peak amplitude, and a 40-50 
days (quasi-)period (Fig. |5]l. To verify that this variability does 
not actually originate in one of the reference stars, we repeated 
the analysis using each of the two reference stars. Those alter- 
nate light curves are very similar to Fig. |5] The variations are 
fully consistent with the 38.92 days period identified by |Kiraga 
& Stepien ( 2007 ) in a much longer photometric timeseries. Our 
photometry demonstrates that Gl 176, which Kiraga & Stepien 



( 2007) find did not significantly vary until JD=2453300, has re- 
mained strongly spotted until the end of our radial velocity mea- 
surements. Our dense sampling also excludes that 38.92 days 
would have been an alias of the true period. We adopt the better 
defined Kiraga & Stepien ( 2007| l value as the rotation period of 
Gl 176. 

Our photometric observations are consistent with the signal 
of a single spot, within the limitations of their incomplete phase 
coverage: the variations are approximately sinusoidal, and their 
~0.2-0.3 phase shift from the corresponding radial velocity sig- 
nal closely matches the difference expected for a spot. The spot 
would cover 2.6% of the stellar surface if completely dark, cor- 
responding to a ~ 0.167?* radius for a circular spot. 

4.2. Variability of the spectroscopic Indices 

The emission reversal in the core of the Ca ii H&K resonant lines 
and in the Ha line results from non-radiative heating of the chro- 
mosphere, which is magnetically coupled to the photospheric 
spots and plages. We measured in the HARPS spectra the spec- 
tral indices defined by Bonfils et al .1 ( 120071 to probe these chro- 
mospheric spectral features, and examine here their variability. 

The power spectra for both the H-hK and Ho- indices have 
clear peaks near 40 days (Fig. |6] lower right panel). Within the 
combined uncertainties, these peaks are consistent with both the 
photometric period and the longer radial velocity period. The 
phasing of the chromospheric index and the photometry is such 
that lower photometric flux matches higher Can emission, as 
expected if active chromospheric regions hover above dark pho- 
tospheric spots. 



Though certainly not as clearly as for Gl 674 (Bonfils et al. 



2007 ), a plot of the (apparent) radial- velocity (after subtraction 
of the 8.8-day planet) against the Hh-K spectral index (Fig.|6] up- 
per panel) similarly suggests the loop pattern which is expected 
for a spot (IBonfils et al. 12007 ): a spot produces maximal velocity 
offsets when it is on either edge of the star, where geometric pro- 
jection reduces the apparent area of its associated chromospheric 
emission to an intermediate value; it produces no velocity offset 
when it crosses the sub-observer meridian, with a maximal pro- 
jected area for a front-facing crossing and a minimal (null for a 
non-polar spot) projected area for a back facing crossing. The ra- 
dial velocity offset therefore cancels for both the minimum and 
the maximum chromospheric emission, and is maximal for in- 
termediate chromospheric emission levels. The pattern here is 
definitely noisier than observed on Gl 674, suggesting that the 
spot pattern may evolve on a time scale of the order of our ob- 
serving period. 
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Fig. 6. Upper panel: Differential radial velocity of Gl 176, corrected for the signature of the 8.8 days planet in our 2-planet fit, as 
a function of the Ha (red filled circles) and Ca ii H&K (green filled squares) spectral indices defined in the text for the 2007/2008 
observing season. Bottom right panels: the Ca ii H+K and Ha indexes phased to the longer period of the 2-planet model. Bottom 
left panels: Power Density spectra of the spectroscopic indices. A clear power excess peaks at 40 days (vertical dashed lines) 



4.3. Planets vs. activity 

In ^ we showed that our 57 radial-velocity measurements of 
Gl 176 are well described by two Keplerian signals. Section]?] 
however demonstrates that the rotation period of Gl 176 coin- 
cides with the longer of these two Keplerian periods. The stellar 
flux and the Can H-hK emission vary with that period, with a 
phase relative to the velocity variations consistent with a mag- 
netic spot on the stellar surface. As a consequence, some, and 
probably all, of the 40-day radial-velocity signal must originate 
in the spot. Planet-induced activity through magnetic coupling 
(e.g. Shkolnik et al.|2005 1 would in principle be an alternative 
explanation of the correlation, but has never been observed for 
such a long-period planet. The inner planet in addition is not 
very much less massive than the hypothetical 40-day planet. One 
would, at least naively, expect its position in the inner magneto- 
sphere of Gl 176 to more than make up for its lower mass. The 
8.8-day period however is only seen in the radial velocity signal, 
and it has no photometric or chromospheric counterpart. 



5. Discussion and conclusions 

The most important result of the above analysis is that a 
M sin(i)=8.4 MEarih planet orbits Gl 176 in a ~8. 8-day orbit. 
Variability identifies the stellar rotation period as 38.92 days, and 
the 8.8-day period therefore cannot reflect rotation modulation. 
The short period signal, in spite of its similar amplitude, also has 



no counterpart in either photometry or chromospheric emission, 
further excluding a signal caused by magnetic activity. 

Like Gl 674 ( |Bonfils et al.|2007| l, Gl 176 demonstrates that 
single planets can be identified around moderately active M- 
dwarfs, at the cost of doubling or tripling the number of measure- 
ments over a magnetically quiet M-dwarf. Since the Keplerian 
model does not reflect a physical reality for the 40-day period, its 
residuals must be interpreted with caution. They are well above 
the measurement errors (x^ - 5.86 per degree of freedom), and 
could in principle reflect addition planet(s) in the system. More 
likely, much of these residuals stem from long term evolution of 
the spot pattern of Gl 176. Many additional radial velocity mea- 
surements would be needed to firmly identify additional planets 
amongst this spot-evolution noise. That cost may, in practice if 
not in theory, effectively impede the detection of multi-planet 
systems around moderately active stars. It may therefore not be 
fully by coincidence that Gl 674 and Gl 176 are both the only 
M-dwarf planetary systems with no indications of further plan- 
ets and the two most active M-dwarfs with known planets. 

At 0.066 AU from its parent star, the thermal equilibrium 
temperature of Gl 176 b is -450 K. Its 8.4 MEanh M sin(i) might 
be sufficient for accretion of a significant gas envelope to have 
occured, in particular in case the inclination would turn out to 
be non-trivial, but the rocky core most likely dominate the total 
mass (e.g. Seager et al.|2007 Valencia et al.|200 7 ). 

With a mass of onl y M sin(i)=8.4 MEar th, Gl 176b adds to the 
growing evidence (e.g. |Bonfils et al.|2007[ ) that super-Earths are 
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common around very low mass stars: 6 of the 20 known planets 'Gl 176' on page 4 

with M sin(i) < 0. 1 Mjup orbit an M dwarf, in contrast to just 3 'Gl 176' on page 4 

of the ~250 known Jupiter-mass planets. 'Gl 674' on page 4 

'Gl 674' on page 4 
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Table 3. Radial-velocity measurements and error bars for 
Gl 176. All values are relative to the solar system barycenter, 
and corrected from the small perspective acceleration using the 
Hipparcos parallax and proper motion. Only available electroni- 
cally. 
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